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Abstract 

We show that the simplest assumptions for the dynamics of particle produc- 
tion allow us to understand the fluxes of hadrons and photons at mountain 
altitudes as well as the structure of individual events. The analysis requires 
a heavy nuclear component of primary cosmic rays above the "knee" in the 

spectrum with average mass number < A >= 7.3 it 0.9. 
PACS number(s): 13.85.Tp, 96.40.De 
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I. INTRODUCTION 



The energy spectrum and chemical composition of primary cosmic rays have been de- 
termined from direct observation above the earth's atmosphere using balloons and space- 
craft IQ. The technique is limited by the small size of the detectors and the short exposure 
times. As a result of the steep energy spectrum no direct observations are available above a 
primary energy of roughly 10^^ eV. In the interesting region of the "knee" in the spectrum 
and above information on composition has to be inferred from indirect measurements of air 
showers at sea level or mountain altitudes 0, by using large area detectors for long periods 
of time. 

In this paper we infer the composition of the cosmic rays from measurements at moun- 
tain altitudes of the hadronic and electromagnetic component of the air cascades initiated 
at the top of the atmosphere. A connection between the nature of the primary particle 
and air shower observations requires the detailed understanding of particle interactions at 
very high energies and forward scattering angles where no information is available from 
accelerator-based experiments. The basic problem is that one is faced with the impossibil- 
ity of deducing two unknowns, the composition and the dynamics of particle interactions, 
from a single measurement. We nevertheless pursue this challenge because we are confident 
that we understand particle interactions with sufficient accuracy to meaningfully approach 
this problem. We have indeed formulated a model which is based on the most straight- 
forward assumptions and which respects the spirit of quantum chromodynamics More 
importantly, it describes in quantitative detail single events, i.e. shower cores in their early 
stage of development, observed in emulsion chamber experiments at mountain altitudes 
Here we will show that this model describes the observed hadronic and electromagnetic 
spectrum at mountain altitudes, provided the mass number of the primary cosmic rays is 
< A >= 7.3 ± 0.9. This is consistent with the result obtained by other indirect means. 

Our model of particle production 0] is guided by the features of QCD-inspired models: 
approximate Feynman scaling in the fragmentation region and an inelasticity slowly varying 
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with energy. The rapidity density of secondary charged pions is parametrized as 

dN dN (1-xY 
dy dx X 

where y is the rapidity of the secondaries and the Feynman variable x is given by the ra- 
tio of the energy E of the secondary particle to the incident energy Eq. With a = 0.12 
and n = 2.6 (n = 3 is expected on the basis of counting rules) the overall features of the 
hadronic component of single events detected in emulsion chambers were quantitatively re- 
produced. For illustration, we present in Fig. the hadronic integral spectrum of two events 
detected by the Brazil- Japan Collaboration at Mt. Chacaltaya, Bolivia (atmospheric depth 
540 g/cm^) Pfl, which are successfully described by our model [0. In the present paper 
we use the same model to calculate both the hadronic and electromagnetic integral spectra 
of atmospheric showers detected in large emulsion chamber experiments. The calculation is 
performed by solving the cosmic-ray diffusion equations using the rapidity distribution for 
particle production given by Eq. ([1|). Starting with the measured all-particle primary spec- 
trum at the top of the atmosphere, we propagate the particle showers down to the mountain 
altitude detection levels of the various experiments and investigate our results as a function 
of the assumed average composition of the primary cosmic radiation. 



II. HADRONIC AND ELECTROMAGNETIC SHOWERS IN THE ATMOSPHERE 

The flux of cosmic ray nucleons at the top of atmosphere (depth t = 0) is parametrized 
by a power-law spectrum 

F„(E,t = 0) = iVoE-(^+i) . (2) 

At this point no secondaries have been produced, hence the boundary condition for the 
pionic component of the shower: F^^lEji = 0) = 0. The hadronic flux Fh{E,t) can be 
calculated for any depth t = z in terms of the interaction mean free path (MFP) of nucleons 
(n) and pions (vr), \i{E) with i = n,n. The result is 
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FhiE,z) = Fn{E,z)+F^{E,z) 



(3) 



with 



z) = No e-"^"(^) , 



(4) 




(5) 



The dependence of the functions Hi{E) and gi{j,E) on the rapidity distribution, Eq. (H), 
and on the MFP Xi{E), is described in the Appendix. 

The electromagnetic component of the shower is initiated by gamma rays from the decay 
of the neutral pion, 7t° 27. With equal multiplicity of vr^, tt^ and 7r° secondaries, the 
number of neutral pions is half the number of the charged pions which is given by Eq. (^. 
The gamma-ray distribution is given by two-body decay 



A gamma ray with energy E^ at depth z contributes to the electromagnetic cascade 
F^{E,t) with 



Here (e + 'y){E^,E,t — z) represents the photons and e+e pairs in the cascade produced 



formalism [|r^. The result is of the form of Eq. (|^) with (e + '-f){E^, E,t — z) given by the 
eigenvalues of the electromagnetic cascade equations; see Appendix. 

We are now ready to compute the integral energy spectrum for hadrons and photons 
{i = /i,7). 






X (e + -f){E^,E,t- z) . 



(7) 



by the photon with energy E^. We compute it in approximation A using the operator 




(8) 



which can be confronted with experimental results. 



4 



III. INELASTIC CROSS SECTION FOR HADRON-AIR INTERACTIONS 



Before proceeding with the analysis, it is necessary to describe the energy dependence 
of the MFP which is inversely proportional to the inelastic cross section for hadron-air 
interactions, i.e. 

24, 100(g/cm2) 



HE) 



a, 



[in mb) 



(9) 



We calculated the inelastic cross section using the event generator SIBYLL |[T^ . The result 
can be parametrized by 



i-air 

^ 1 n ^7, 



1 + 6, In^ 



E 
E'o 



(10) 



with Eq = 200 GeV. Our results with s„ = 284.5 mb and 6„ = 3.852 x 10 ^ for proton-air 
and = 211.0 mb and fe^r = 5.827 x 10^'^ for pion-air scattering are shown in Fig. S. 



IV. PRIMARY COMPOSITION AND SHOWER ENERGY SPECTRA 

Having constructed an explicit model of particle interactions which successfully describes 
individual events, see Fig. 0, we can compute the flux of hadrons and photons at mountain 
altitude as a function of the primary cosmic ray flux. For the primary spectrum we use a 
parametrization |l^ which is accurate in the energy region between 300 and 10^ TeV/particle 
relevant to our calculation. It accurately extrapolates to lower energy measurements ob- 



tained with the Proton satellite and the JACEE balloon flights [p!3| , p!4| . The all-particle 
differential energy spectrum, is given by 



Faii{E,t = 0) = (4.55 ±0.45) x 10 



-11 



(11) 



103.67_ 

in units of (m^ sec sr TeV/particle) with 7 = 1.62±0.12 below and 7 = 2.02±0.05 above 
TeV/particle. Eq. (^) describes the change in the slope of the spectrum at the energy 
region known as "the knee". We parametrize our ignorance of the chemical composition of 
the primary flux in terms of a single parameter < A >, the average mass number of the 



primary nuclei. Heavy primaries are included in our formalism using superposition in 
Eq. (^). The projectile nucleus of energy Eq is considered to be the superposition of A 
nucleons interacting independently, each having energy Eq/A. Although a simplification, 
this assumption is quite acceptable, as long as the nuclei fragment relatively rapidly. 

We first calculate the integral energy spectra of electromagnetic showers, Eq. (|), at the 
detection level of Mt. Chacaltaya, using the extremes values for < A > corresponding to 
pure proton {A = 1) and to pure iron {A = 56). The results are shown in Fig. ^, by the 
dashed curves. The predictions bracket the experimental data [l^. That a pure proton 
spectrum cannot describe these results is not totally surprizing It is well known that 
the relatively low rate of detected gamma-ray families (and also of halo families) cannot be 
understood in models with approximate Feynman scaling unless heavy primaries contribute 
to the cosmic ray flux. We subsequently determine, by chi-square minimization, the average 
mass number that best describes the data. We obtain < A >= 7.3 ±0.9 (solid line in Fig. ^. 

Having fixed all parameters, we can confront the model with any other observations. We 
find that it successfully describes both the hadronic (Fig. ^) and electromagnetic (Fig. ^) 
components of the atmospheric showers detected in large emulsion chambers at Mt. Fuji 
in Japan (atmospheric depth 650 g/cm^), and at Mt. Kanbala |]TB| in China (520 g/cm^). 



V. CONCLUSIONS 

We conclude that with the simplest assumptions for the production of secondaries based 
on approximate scaling in the fragmentation region, it is possible to explain a broad set of 
experimental data on very high energy cosmic rays in the atmosphere, namely the lateral 
spread and the integral spectra of superfamilies (as in Ref. and Fig. |l|), and the energy 
spectra of hadronic and electromagnetic showers detected in large emulsion chamber experi- 
ments (as in Figs, ^and This scenario requires a primary composition with average mass 
number 7.3 ± 0.9. We investigated the sensitivity of this quantity to different parametriza- 
tions of the all-particle spectrum [0,^] and the best fit invariably yields < A >~ 7. Our 
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result is also consistent with underground muon measurements which yield an A-value 
of 10 ± 4 in the 1 to 1,000 TeV energy range. 



It has been noted elsewhere |^ that it is difficult to establish whether one must adopt 
a heavy primary composition along with a model of particle production based on scaling 
or, alternatively, a proton dominant composition along with a strong violation of Feynman 
scaling. It should be noted however that in our analysis the particle interaction model was 
determined on the basis of an independent study of single events initiated by protons deep 
in the atmosphere. The a posteriori analysis of the hadron and photon spectra at mountain 
altitude presented here, required the introduction of heavy primaries. 
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APPENDIX: DEFINITIONS IN THE ANALYTICAL SOLUTIONS 

Complete definition of the hadronic flux components presented in Eqs. (H) and re- 
quires the following functions (for z = ra, tt): 



1 / "i 



K(EI 



(Al) 



/•I 1 dN 



with 
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/■I dN 



W7) = / f{(Ti)(T]dai , 

JO 

\ /■! 1 



f{ai)a]d(ri 



\\i{E/ai)l Jo \{E/a, 
The elasticity distribution is assumed to be 

/(a,) = (l + /?) , (A4) 

where (3 fulfills a consistency relation between average elasticity (cr) and average inelasticity 
{K), so that (cr) + [K) = 1 (energy conservation). The eigenvalues for the electromagnetic 
cascade equations are 

1 Bjs) 
^o(Ai(s)-A2(s)) 

X (e^iWt/Xo_gA,(s)VXo^ ^ (A5) 



+ e^^W*/^«) , (A6) 



where Hi{s), H2{s), Ai(s), A2(s) and Xq are parameters with standardized definitions in cas- 
cade theory ||10|- Subsequently (e + 'y){E^,E,t — z) in Eq. (|^) should be replaced by 
(n^(s, t — z) +V^{s,t — z)) with s evaluated at the pole 5 = 7. 
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FIGURES 

FIG. 1. Integral energy spectra of hadronic superfamily events detected at Mt. Chacaltaya [5,6]. 
Ursa Maior event (0) and CentauroVII data (A), are compared to the calculation of Rcf. [3] (solid 
line), using the x-distribution of Eq. (1). For illustrative purposes, the data of CentauroVII have 
been shifted by a factor 10. 

FIG. 2. Inelastic cross sections for p— air (Q) and tt— air (□) scattering computed using 
SIBYLL [12] and parameterized by Eq. (10). 

FIG. 3. Integral energy spectra of electromagnetic showers (0), detected at Mt. Chacaltaya [15], 
compared to the calculation using the x-distribution of Eq. (1). Dashed lines: < A >= 1 (proton) 
and 56 (iron); solid line: < A >= 7.3 it 0.9, dotted lines: calculated from uncertainties in A and 7. 

FIG. 4. Integral energy spectra of showers detected at Mt. Fuji [17] (Q) and Mt. Kan- 
bala [18] (A), compared to the analytical calculation using the a;-distribution of Eq. (1), with 
< A >= 7.3 lb 0.9 (solid line), (a) Hadron induced showers; (b) Electromagnetic showers. The data 
of Mt. Kanbala have been shifted by a factor 100. Dotted lines are calculated from uncertainties 
in A and 7. 



11 




E'^ (TeV) 



Fig. 1 



1000 




Fig. 2 




Fig. 4a 




Fig. 4b 



